Chapter 5 



Electromagnetic Waves in Plasmas 



5.1 General Treatment of Linear Waves in Anisotropic 
Medium 

Start with general approach to waves in a hnear Medium: Maxwell: 

1 (9E (9B 
VAB = „-+^^ ; VAE=-^ (5.1) 

we keep all the medium's response explicit in j. Plasma is (infinite and) uniform so we Fourier 
analyze in space and time. That is we seek a solution in which all variables go like 

expi(k.x — Lot) [real part of] (5.2) 

It is really the linearised equations which we treat this way; if there is some equilibrium field 
OK but the equations above mean implicitly the perturbations B, E, j, etc. 

Fourier analyzed: 

ik A B = /ij + -^E ; ik A E = icuB (5.3) 
Ehminate B by taking kA second eq. and u; x 1st 

ik A (k A E) = ujfXoj (5.4) 



So 

kA (kAE) + ^E + ia;yuJ = (5.5) 



a;2 



C2 



Now, in order to get further we must have some relationship between j and E(k, cj). This 
will have to come from solving the plasma equations but for now we can just write the most 
general linear relationship j and E as 

j = 0-.E (5.6) 
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cr is the 'conductivity tensor'. Think of this equation as a matrix e.g.: 

(5.7) 





This is a general form of Ohm's Law. Of course if the plasma (medium) is isotropic (same 
in all directions) all off-diagonal a' s are zero and one gets j = crE. 

Thus 



2 



k(k.E) - /c^E + ^E + iuJUoCr.'E = (5.8) 

Recall that in elementary E&M, dielectric media arc discussed in terms of a dielectric con- 
stant e and a "polarization" of the medium, P, caused by modification of atoms. Then 

e^E = ^ - and V.D = pext (5.9) 

Displa<;ement Polarization externalcharge 

and one writes 

P= ^ e„E (5.10) 
susceptibility 

Our case is completely analogous, except we have chosen to express the response of the 
medium in terms of current density, j, rather than "polarization" P For such a dielectric 
medium. Ampere's law would be written: 

1 dH d 

— VAB=jext+^ = T^eeoE, if jext = , (5.11) 
Ho at ot 

where the dielectric constant would be e = 1 -|- 

Thus, the explicit polarization current can be expressed in the form of an equivalent dielectric 
expression if 

dE ^ dE d ^ 
, + eo- = a.E + e-^-eoe.E (5.12) 

or 

e = 1 + — ^ (5.13) 
—icueo 

Notice the dielectric constant is a tensor because of anisotropy. The last two terms come 
from the RHS of Ampere's law: 

j + |(eoE) . (5.14) 

If we were thinking in terms of a dielectric medium with no explicit currents, only implicit (in 
e) we would write this ^(eepE); e the dielectric constant. Our medium is possibly anisotropic 
so we need ^(eoe.E) dielectric tensor. The obvious thing is therefore to define 

e = 1 + ^— o- = 1 + (5.15) 

—lU!€o CO 
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Then 



k(k.E) - rE + — e.E = 



and we may regard e(k, c<j) as the dielectric tensor. 
Write the equation as a tensor multiplying E: 



D.E = 



with 



UJ 



D = {kk - ri + -^e} 



Again this is a matrix equation i.e. 3 simultaneous homogeneous eqs. for E. 



D. 



yx 



( Ex 

Ey 



In order to have a non-zero E solution we must have 

det I D 1= 0. 



(5.16) 



(5.17) 

(5.18) 

(5.19) 



(5.20) 



This will give us an equation relating k and cu, which tells us about the possible wavelengths 
and frequencies of waves in our plasma. 

5.1.1 Simple Case. Isotropic Medium 



cr — a 1 

e = e 1 

Take k in z direction then write out the Dispersion tensor D. 



D 
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kk / \ 



kk 



A:^ 

A;2 



+ 



4e 



\ 

V ^0 ie) 





73-' 



(5.21) 
(5.22) 



+ 







Take determinant: 

det I D 1= [-1 
Two possible types of solution to this dispersion relation: 





























e = 




J C2 



(5.23) 



(5.24) 
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(A) 



, 2 ^ 



= 0. 













E,. 












Ey 










\ 





E, = 0. 



(5.25) 
(5.26) 



Electric field is transverse (E.k — 0) 
Phase velocity of the wave is 



k 



c 



This is just like a regular EM wave traveling in a medium with refractive index 

kc 



N 



UJ 



(B) 



UJ 





^ E^xx 













Eyy 







I 









-e = i.e. e = 



Er,. 



Ey 

E. 



Ex — Ey — 0. 



(5.27) 

(5.28) 

(5.29) 
(5.30) 



Electric Field is Longitudinal (E A k = 0) E || k. 

This has no obvious counterpart in optics etc. because e is not usually zero. In plasmas 
e = is a relevant solution. Plasmas can support longitudinal waves. 



5.1.2 General Case (k in z-direction) 



D 



^zx 



-N^ + e 

^zy 



yy 



^xz 

^yz 

^zz 



k^c 



(5.31) 



When we take determinant we shall get a quadratic in N"^ (for given uj) provided e is not 
explicitly dependent on k. So for any uj there are two values of A^^. Two 'modes'. The 
polarization E of these modes will be in general partly longitudinal and partly transverse. 
The point: separation into distinct longitudinal and transverse modes is not possible in 
anisotropic media (e.g. plasma with Bo). 

All we have said applies to general linear medium (crystal, glass, dielectric, plasma). Now 
we have to get the correct expression for cr and hence e by analysis of the plasma (fluid) 
equations. 
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5.2 High Frequency Plasma Conductivity 



We want, now, to calculate the current for given (Fourier) electric field E(k, cu), to get the 
conductivity, cr. It won't be the same as the DC conductivity which we calculated before 
(for collisions) because the inertia of the species will be important. In fact, provided 

u > Pet (5.32) 

we can ignore collisions altogether. Do this for simplicity, although this approach can be 
generalized. 

Also, under many circumstances we can ignore the pressure force — Vp. In general will 
be true if ^ 3> fte,i We take the plasma equilibrium to be at rest: Vq = 0. This gives a 
manageable problem with wide applicability. 

Approximations: 

CoUisionless = 

'Cold Plasma' Vp = {e.g. T ~ 0) (5.33) 

Stationary Equil Vo — 

5.2.1 Zero B-field case 

To start with take Bq = 0: Plasma isotropic Momentum equation (for electrons first) 

(9v 



mn 



at + 



nqE (5.34) 



Notice the characteristic of the cold plasma approx. that we can cancel n from this equation 
and on linearizing get essentially the single particle equation. 

m-^ = gE (Drop the 1 suffix now). (5.35) 
This can be solved for given uj as 

V = — ^E (5.36) 

—turn 

and the current (due to this species, electrons) is 

2 

j = ngv = ^^E (5.37) 
—tcum 

So the conductivity is 

(7 = i^ (5.38) 
cum 

Hence dielectric constant is 

6 = l + ^a = l-f^H = l + X (5.39) 
ooeo \meoJ OJ 



100 



Longitudinal Waves (Bq = 0) 

Dispersion relation we know is 

, = 0=1- f^U (5.40) 

\meoJ OJ^ 

[Strictly, the e we want here is the total e including both electron and ion contributions to 
the conductivity. But 

^ ~ ^ (for z = 1) (5.41) 

(Tj me 

SO to a first approximation, ignore ion motions.] 
Solution 

= f (5.42) 
\meeo J 

In this approx. longitudinal oscillations of the electron fluid have a single unique frequency: 

This is called the 'Plasma Frequency ' (more properly ujpe the 'electron' plasma frequency) . 
If we allow for ion motions we get an ion conductivity 



urrii 

and hence 



(5.44) 



etot = IH ((Te + (Ti) = 1 \ — (5.45) 

ujEo yeorrie eorrii I uj^ 



where 



\eomi 



is the 'Ion Plasma Frequency'. 



Simple Derivation of Plasma Oscillations 

Take ions stationary; perturb a slab of plasma by shifting electrons a distance x. Charge 
built up is rieqx per unit area. Hence electric field generated 
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Figure 5.1: Slab derivation of plasma oscillations 



Equation of motion of electrons 

dv n^qlx 

me-7- = ; (5.48) 

at Eg 

Simple harmonic oscillator with frequency 

^ Plasma Frequency. (5.50) 

eome) 

The Characteristic Frequency of Longitudinal Oscillations in a plasma. Notice 

1. a; = u;p for all k in this approx. 

2. Phase velocity | can have any value. 

3. Group velocity of wave, which is the velocity at which information/energy travel is 

(5.51) 

In a way, these oscillations can hardly be thought of as a 'proper' wave because they do 
not transport energy or information. (In Cold Plasma Limit). [Nevertheless they do emerge 
from the wave analysis and with less restrictive approxs do have finite Vg] 

Transverse Waves (So = 0) 
Dispersion relation: 

- A;^ + — e = (5.52) 
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or 



7V2 = 



1 - + ul) lu" 



pe 



(5.53) 




Figure 5.2: Unmagnetized plasma transverse wave. 




Figure 5.3: Alternative dispersion plot. 



Alternative expression: 



-A;^ + '^-4 = (5.54) 
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which imphes 



2 2 I ; 2 2 

uj = a;„ + fee 



(5.55) 
(5.56) 



5.2.2 Meaning of Negative N'^: Cut Off 

When iV^ < (for uj < LVp) this means N is pure imaginary and hence so is A; for real uo. 
Thus the wave we have found goes hke 



exp{± \ k \ X — icut} 



(5.57) 



i.e. its space dependence is exponential not oscillatory. Such a wave is said to be 'Evanescent' 
or 'Cut Off'. It does not truly propagate through the medium but just damps exponentially. 

Example: 




Figure 5.4: Wave behaviour at cut-off. 



A wave incident on a plasma with cUp > uP' is simply reflected, no energy is transmitted 
through the plasma. 



5.3 Cold Plasma Waves (Magnetized Plasma) 

Objective: calculate e,D,k(a;), using known plasma equations. 
Approximation: Ignore thermal motion of particles. 

Applicability: Most situations where (1) plasma pressure and (2) absorption are negligible. 
Generally requires wave phase velocity ^ fthermai- 
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5.3.1 Derivation of Dispersion Relation 

Can "derive" the cold plasma approx from fluid plasma equations. Simpler just to say that 
all particles (of a specific species) just move together obeying Newton's 2nd law: 



(9v 

m— = qCE + V A B) 

ot 



(5.58) 



Take the background plasma to have Eq = 0, B = Bq and zero velocity. Then all motion 
is due to the wave and also the wave's magnetic field can be ignored provided the particle 
speed stays small. (This is a linearization). 



m— =g(E + vABo), 



(5.59) 



where v, E oc expi(k.x — ut) arc wave quantities. 
at 



Substitute ^ — >• —iu and write out equations. Choose axes such that Bq = -Bo(0,0, 1). 



— lumvx 
—iujmv,, 



q{Ey - v^Bq) 



(5.60) 



Solve for v in terms of E. 



iujmVz — qEz 

q (iujEx — ^Ey 



Vcr. = 



m 



00" 



9? 



q / VtEx + iuoEy 

m\ uj"^ — Q? 
q i 



(5.61) 



-E, 



where Vl = ^ is the gyrofrequency but its sign is that of the charge on the particle species 
under consideration. 

Since the current is j = qvn — <t.E we can identify the conductivity tensor for the species 
(j) as: 

















rrij aj2— n| 



















(5.62) 



The total conductivity, due to all species, is the sum of the conductivities for each 



(5.63) 



So 



'yy 



qfrij lUJ 
f rrij up' — 



E 



(5.64) 
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a 



xy 



—a. 



yx 



-E 



, rrii up' — 



Susceptibility x — zr—^^- 
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-iD 





^yx ^yy 
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S 





e,, _ 










p 



where 



"^^xy — 



-yy 



-le. 



yx 



_ _ ^ ^Vd 



1-E 



uo. 



UJ 



and 



P3 



eoTTlj 



is the "plasma frequency" for that species. 
S & D stand for "Sum" and "Difference" : 

S=l{R + L) D = \{R-L) 

where R L stand for "Right-hand" and "Left-hand" and are: 



R^l-Y. 



UJ. 



UJ (uj + flj 



1-E 



uj: 



UJ (uj — 



(5.65) 
(5.66) 



(5.67) 



(5.68) 
(5.69) 
(5.70) 

(5.71) 



(5.72) 



(5.73) 



The R Sz L terms arise in a derivation based on expressing the field in terms of rotating 
polarizations (right & left) rather than the direct Cartesian approach. 

We now have the dielectric tensor from which to obtain the dispersion relation and solve 
it to get k(a;) and the polarization. Notice, first, that e is indeed independent of k so the 
dispersion relation (for given cu) is a quadratic in N"^ (or k"^). 

Choose convenient axes such that ky = Ny — 0. Let 9 be angle between k and Bq so that 

N, = N cose, N^^Nsine. (5.74) 

Then 



D = 



-iV^cos^^-h^ -iD sine cos 9 

+iD -N^ + S 

N'^sm9cos9 -N^sm^9 + P 



(5.75) 
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and 
where 



Solutions are 



WD \\^ AN^ - BN^ + C (5.76) 

A = S sin^ e + P cos'^ e (5.77) 

B = RLsin^e + PS{l + cos'^e) (5.78) 

C = PRL (5.79) 



- (5.80) 



where the discriminant, F, is given by 

= {RL - PSf sin^ 9 + AP'^D'^ cos^ 9 (5.81) 

after some algebra. This is often, for historical reasons, written in the equivalent form (called 
the Appleton-Hartree dispersion relation) 

The quantity is generally +ve, so is real =^ "propagating" or "evanescent" no wave 
absorption for cold plasma. 

Solution can also be written 

2^ P (N^ - R) (N^ - L) 

This compact form makes it easy to identify the dispersion relation at 6* = & | i.e. parallel 
and perpendicular propagation tan^ = 0, oo. 

Parallel: P = Q , = R N'^ = L 

Perp: N'^ = ^ N'^ = P . 

Example: Right-hand wave 

N"^ — R. (Single Ion Species). 

2 2 

, "^"^ „ - , "^^^1^ I, (5.84) 

This has a wave resonance N"^ — > oo at a; = |r2e|, only. Right-hand wave also has a cutoff at 
i? = 0, whose solution proves to be 



I I I ^^'j I 

CJ = UJR= h 



I \ 2 

I I ( 



I + + a;pj 



1/2 

(5.85) 
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Since ^ mg this can be approximated as: 

I a 



1+1 + 4- 



uo. 



pe 



This is always above \Q,e 




(5.86) 



Figure 5.5: The form of the dispersion relation for RH wave. 

One can similarly investigate LH wave and perp propagating waves. The resulting wave 
resonances and cut-offs depend only upon 2 properties (for specified ion mass) (1) Density 
ujp^ (2) Magnetic Field ^ |0e|- [Ion values ujpi, \VLi\ are got by ^ factors.] 

These resonances and cutoffs are often plotted on a 2-D plane ^ , ^ (oc n) called the 

C M A Diagram. 

We don't have time for it here. 



5.3.2 Hybrid Resonances Perpendicular Propagation 

"Extraordinary" wave N"^ — ^ 



^2 _ 02) (^2 _ _ ^2 (^2 _ J^2) _ (^2 _ J^2) 



Resonance is where denominator = 0. Solve the quadratic in u"^ and one gets 



00, 



pe 



02 



(5.87) 



(5.88) 
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Neglecting terms of order ^ (e.g. ^) one gets solutions 



ojjjff — cUpg + Upper Hybrid Resonance. 



UJ 



LH 



e pi 



02 + 



2 

pe 



Lower Hybrid Resonance. 



At very high density, ujp^ 3> fl 



geometric mean of cyclotron frequencies. 
At very low density, ujp^ <C Og 

ion plasma frequency 

Usually in tokamaks cu^g ~ Ql- Intermediate. 
Summciry Graph > Up) 



(5.89) 
(5.90) 

(5.91) 
(5.92) 




Figure 5.6: Summary of magnetized dispersion relation 

Cut-off's arc where = 0. 
Resonances are where N'^ — ^ oo. 

Intermediate angles of propagation have refractive indices between the 6 — 0, ^ lines, in the 
shaded areas. 
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5.3.3 Whistlers 



(Ref. R.A. Helliwell, "Whistlers & Related Ionospheric Phenomena," Stanford UP 1965.) 
For N"^ ^ 1 the right hand wave can be written 

,2 



Group velocity is 



Then since 



we have 



(Ncu) 



Thus 



Group Delay is 



—OJ. 



pe 



CO {CU — \ fle\) ' 

duo f dk\ ^ 



{N = kc/uj) 



d /Nuj\ 
duj \ c J 



n -1 



duj ( 






UJp/2 




— Uj)'^ UJ2 




p 1 1 /2 




— a;j 2 a;2 




c 



1 ) 
2 



-{(l^el -'^)+a;} 



c2i\n,\ -^Yuj2 



Up\Q,e\ 



— oc — 3- oc ^ 



(5.93) 

(5.94) 
(5.95) 



(5.96) 

(5.97) 
(5.98) 




Figure 5.7: Whistler delay plot 
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Plot with — as a;-axis. 
Resulting form explains downward whistle. 
Lightning strike ~ 5-function excites all frequencies. 
Lower ones arrive later. 



Examples of actual whistler sounds can be obtained from http://www-istp.gsfc.nasa 



gov/istp/polar/polar_pwi_sounds . html 



5.4 Thermal Effects on Plasma Waves 

The cold plasma approx is only good for high frequency, N"^ ~ 1 waves. If u is low or A^^ ^ 1 
one may have to consider thermal effects. From the fluid viewpoint, this means pressure. 
Write down the momentum equation. (We shall go back to Bq = 0) linearized 

mn—— = nqEi — Vpi ; (5.99) 
remember these are the perturbations: 

P = Po+Pi ■ (5.100) 

Fourier Analyse (drop I's) 

mn{—iuj)\^ = ngE — ikp (5.101) 
The key question: how to relate p to v 
Answer: Equation of state + Continuity 
State 

pn~'' = const. =^ {pq + pi) {uq + ni)~^ = PqUq'^ (5.102) 
Use Taylor Expansion 



(Po + Pi) {no + ni) ~ PoUq ^ 

Hence 



, , Pi ni 

IH 7 — 

Po no 



Continuity 



Linearise: 



Fourier Transform 



(5.103) 



^ = 7^ (5.104) 

Po no 



— + V. (r2v)=0 (5.105) 
at 

^ + V. (novi) = — + r2oV.v = (5.106) 

— iuui + noik.^i = (5.107) 



111 



I.e. 



k.v 



ni = no- 



UJ 



Combine State & Continuity 



ni rio 

Pi =Po1— =Po7 

no Uo 



k.v 



k.v 



Po7- 



Hence Momentum becomes 



mn {—^oj) V = nqhi k.v 



(5.108) 



(5.109) 



(5.110) 



Notice Transverse waves have k.v = 0; so they are unaffected by pressure. 

Therefore wc need only consider the longitudinal wave. However, for consistency let us 
proceed as before to get the dielectric tensor etc. 

Choose axes such that k = ke^ then obviously: 



iq 



" cum " 



m —iuj + (ik^^po/mnuj) 



Hence 



mq 
cum 



1 
1 










e = l + 



eoOJ 



1 - ^ 





1 












(5.111) 
(5.112) 

(5.113) 
(5.114) 



mn 

(Taking account only of 1 species, electrons, for now.) 

We have confirmed the previous comment that the transverse waves (E^jEy) are unaffected. 
The longitudinal wave is. Notice that e now depends on k as well as cu. This is called 'spatial 

dispersion'. 

For completeness, note that the dielectric tensor can be expressed in general tensor notation 
as 

,2 / 



LU^ 



1 + kk 



1 - 



1 



- 1 



= 1 



1 + kk- 



(5.115) 



This form shows isotropy with respect to the medium: there is no preferred direction in 
space for the wave vector k. 



112 



But once k is chosen, e is not isotropic. The direction of k becomes a special direction. 
Longitudinal Waves: dispersion relation is 

^zz — (as before) (5.116) 

which is 

1 — = . (5.117) 

mn 



or 

mn 



u^ = ul + e^ (5.118) 



Recall po = UqT = uT; so this is usually written: 

^ul + k^'^ ^ul + k^-fv^ (5.119) 



m 



[The appropriate value of 7 to take is 1 dimensional adiabatic i.e. 7 = 3. This seems plausible 
since the electron motion is 1-d (along k) and may be demonstrated more rigorously by kinetic 
theory.] 

The above dispersion relation is called the Bohm-Gross formula for electron plasma waves. 
Notice the group velocity: 

duo 1 duo'^ 'ykVf 
^9 = -n- = ^.-n- = . * 1 7^ 0. (5.120) 



dk 2u! dk 



and for kvt > ojp this tends to 72^^. In this limit energy travels at the electron thermal 
speed. 

5.4.1 Refractive Index Plot 

Bohm Gross electron plasma waves: 

= 74 fl - ^1 (5.121) 



Transverse electromagnetic waves: 



7eVte V 



Ar2 = _ (5.122) 

These have just the same shape except the electron plasma waves have much larger vertical 
scale: 

On the E-M wave scale, the plasma wave curve is nearly vertical. In the cold plasma it was 
exactly vertical. 

We have relaxed the Cold Plasma approximation. 
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Figure 5.8: Refractive Index Plot. Top plot on the scale of the Bohm-Gross Plasma waves. 
Bottom plot, on the scale of the E-M transverse waves 



5.4.2 Including the ion response 

As an example of the different things which can occur when ions are allowed to move include 
longitudinal ion response: 

up' up 

^ ^ ^jj2 _ k'^Pele ^2 _ k'^Pili (5.123) 

meUe rmrii 

This is now a quadratic equation for cj^ so there are two solutions possible for a given uj. One 
will be in the vicinity of the electron plasma wave solution and the inclusion of uj^^ which is 
^ cUpg will give a small correction. 

Second solution will be where the third term is same magnitude as second (both will be 
^ 1). This will be at low frequency. So we may write the dispersion relation approximately 
as: 



col u;^ 



I.e. 



UJ — 

rriene rrnni 



2 _ k^pai (^li k'^Pele 



=0 (5.124) 



UJ^ = + 
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liPi ^ JePe 
Hi Ue 
liTi + ^eTe 



1 



rrii 



(5.125) 



[In this case the electrons have time to stream through the wave in 1 oscillation so they 
tend to be isothermal: i.e. 7e — 1. What to take for 7^ is less clear, and less important 
because kinetic theory shows that these waves we have just found are strongly damped 
unless Ti < Te] 

These are 'ion-acoustic' or 'ion-sound' waves 



Cg is the sound speed 



— — c 



rrii 



rrii 



(5.126) 



(5.127) 



Approximately non-dispersive waves with phase velocity Cg- 



5.5 Electrostatic Approximation for (Plasma) Waves 

The dispersion relation is written generally as 

N A (N A E) + e.E = N(N.E) - A^^E + e.E = (5.128) 

Consider E to be expressible as longitudinal and transverse components E^, E^ such that 
NAE^ = 0, N.Et = 0. Then the dispersion relation can be written 

N (N.E^) - TV^ + Ej) + e. (E^ + E^) = -N'^Y^t + e.Ej + e.Et = (5.129) 

or 

(Ar2 _ .Et = e.Ee (5.130) 

Now the electric field can always be written as the sum of a curl-free component plus a 
divergenceless component, e.g. conventionally 

E= -V(/) + ,_A^ (5.131) 

Curl-free Divergence- free 

Electrostatic Electromagnetic 

and these may be termed electrostatic and electromagnetic parts of the field. 

For a plane wave, these two parts are clearly the same as the longitudinal and transverse 
parts because 

— V0 = — ik0 is longitudinal (5.132) 
and if V.A — (because V.A — (w.l.o.g.)) then k.A = so A is transverse. 
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'Electrostatic ' waves are those that are describable by the electrostatic part of the electric 
field, which is the longitudinal part: \Ee\ ^ \Et\. 

If we simply say = then the dispersion relation becomes e.E^ = 0. This is not the most 
general dispersion relation for electrostatic waves. It is too restrictive. In general, there is 
a more significant way in which to get solutions where \Ei\ ^ \Et\. It is for iV^ to be very 
large compared to all the components of e : A^^ e ||. 

If this is the case, then the dispersion relation is approximately 

N^Et^e.Ei ; (5.133) 

Et is small but not zero. 

We can then annihilate the E^ term by taking the N component of this equation; leaving 

N.e.E^ = (N.e.N) Ee ^ : k.e.k = . (5.134) 



When the medium is isotropic there is no relevant difference between the electrostatic dis- 
persion relation: 

N.e.N = (5.135) 

and the purely longitudinal case e.N = 0. If we choose axes such that N is along z, then the 
medium's isotropy ensures the off-diagonal components of e are zero so N.e.N = requires 
^zz — ^ e.N = 0. However if the medium is not isotropic, then even if 

N.e.N (= A^^e^^) = (5.136) 

there may be off-diagonal terms of e that make 

e.N ^ (5.137) 

In other words, in an anisotropic medium (for example a magnetized plasma) the electrostatic 
approximation can give waves that have non-zero transverse electric field (of order ||e||/iV^ 
times E£) even though the waves are describable in terms of a scalar potential. 

To approach this more directly, from Maxwell's equations, applied to a dielectric medium 
of dielectric tensor e, the electrostatic part of the electric field is derived from the electric 
displacement 

V.D = V. (eoe.E) = p = (no free charges) (5.138) 
So for plane waves = k.D = k.e.E = ik.e.k0. 

The electric displacement, D, is purely transverse (not zero) but the electric field, E then 
gives rise to an electromagnetic field via V A H = dY)/dt. If A^^ ^11 ^ II then this magnetic 
(inductive) component can be considered as a benign passive coupling to the electrostatic 
wave. 

In summary, the electrostatic dispersion relation is k.e.k = 0, or in coordinates where k is 
in the z-direction, Ezz — ^- 
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5.6 Simple Example of MHD Dynamics: Alfven Waves 



Ignore Pressure & Resistance. 



Linearize: 



E + V AB = 

V = Vi, B = Bo + Bi (Bo uniform), j = ji. 
(9V 



P 



dt 



J A Be 



Fourier Transform: 



E + V A Bo = 

p(-zcu)V = j A Bo 
E + V A Bo = 



Eliminate V by taking 5.144 ABo and substituting from 5.145 

E + ^— (j A Bo) A Bo = 



or 



E 



-tLUp 



^ -{(j.Bo)Bo-i?o'j}^ ^° 



So conductivity tensor can be written {z in B direction). 



R2 

-Do 



1 
1 
oo 



where oo implies that E\\ = (because of Ohm's law). Hence Dielectric Tensor 



e = 1 + 



—tueo 



1 + 



eoi?^ 



1 
1 
oo 



Dispersion tensor in general is: 



D 



NN - A^^ + e 



(5.139) 
(5.140) 

(5.141) 
(5.142) 
(5.143) 

(5.144) 
(5.145) 

(5.146) 
(5.147) 



(5.148) 



(5.149) 



(5.150) 



Dispersion Relation taking A^^ = N^, Ny = 



|D| 

















oo 



(5.151) 
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T~~' T 1 — 




C 



Figure 5.9: Compressional Alfven Wave. Works by magnetic pressure (primarily). 



Meaning of oo is that the cofactor must be zero i.e. 



-m + 1 



-N^ + 1 + 







(5.152) 



The I's here come from Maxwell displacement current and are usually negligible {N]_ ^ 1). 
So final waves are 



1. N'^ = =^ Non-dispersive wave with phase and group velocities 



where we call 

= va the 'Alfven Speed' 

Polarization: 

Eii^E.^O, E^^O. Ey^O ^Vy^O T4^0 (T4 = 0) 



(5.153) 



(5.154) 



(5.155) 

Party longitudinal (velocity) wave — > Compression "Compressional Alfven Wave". 

2. 

II eo-B^ a;"' 

Any 00 has unique Wave has unique velocity in || direction: va- 
Polarization 



E.^Ey^O E^^O ^14 = Vyy^O {V,^0) 
Transverse velocity: "Shear Alfven Wave". 

Works by field hne bending (Tension Force) (no compression). 



(5.156) 
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Figure 5.10: Shear Alfven Wave 



5.7 Non-Uniform Plasmas and wave propagation 

Practical plasmas are not infinite & homogeneous. So how does all this plane wave analysis 
apply practically? 

If the spatial variation of the plasma is slow c.f. the wave length of the wave, then coupling 
to other waves will be small (negligible). 



Figure 5.11: Comparison of sudden and gradualy refractive index change. 

For a given uj, slowly varying plasma means N/^ A or kN/^ ^ 1. Locally, the plasma 
appears uniform. 

Even if the coupling is small, so that locally the wave propagates as if in an infinite uniform 
plasma, we still need a way of calculating how the solution propagates from one place to 
the other. This is handled by the 'WKB(J)' or 'eikonal' or 'ray optic' or 'geometric optics' 
approximation. 

WKBJ solution 

Consider the model 1-d wave equation (for field lu) 





(5.157) 



with k now a slowly varying function of x. Seek a solution in the form 

E — exp (i0 (x)) {—iujt implied) 



(5.158) 
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(f) is the wave phase {— kx in uniform plasma). 
Differentiate twice 



dx^ dx^ \ dx 
Substitute into differential equation to obtain 

/ 7 , \ 2 



2 , ..rfV 



Recognize that in uniform plasma ^ = 0. So in slightly non-uniform, 1st approx is to 
ignore this term. 

^ ~ ±k(x) (5.161) 
dx 

Then obtain a second approximation by substituting 



so 



~ e±i^ (5.163) 



dx J dx 
dx \ 2k dx 



i dk\ 

~ ± Ik ± ——] using Taylor expansion. (5.164) 



Integrate: 



Hence E is 



(^^±±1 kdx + ih\{k^^^ (5.165) 



E^e'^ = cxp (^±i J'' kdx^ (5.166) 

This is classic WKBJ solution. Originally studied by Green & Liouville (1837), the Green 
of Green's functions, the Liouville of Sturm Liouville theory. 

Basic idea of this approach: (1) solve the local dispersion relation as if in infinite homogeneous 
plasma, to get k{x), (2) form approximate solution for all space as above. 

Phase of wave varies as integral of kdx. 

In addition, amplitude varies as -\. This is required to make the total energy flow uniform. 

A;2 



5.8 Two Stream Instability 

An example of waves becoming unstable in a non-equilibrium plasma. Analysis is possible 
using Cold Plasma techniques. 

Consider a plasma with two participating cold species but having different average velocities. 
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These are two "streams" . 

Speciesl Species2 

Moving. Stationary. (5.167) 
Speed V 

We can look at them in different inertial frames, e.g. species (stream) 2 stationary or 1 
stationary (or neither). 

We analyse by obtaining the susceptibihty for each species and adding together to get total 
dielectric constant (scalar 1-d if unmagnetized) . 

In a frame of reference in which it is stationary, a stream j has the (Cold Plasma) suscepti- 
bility 

X. = ^ . (5.168) 

UJ 

If the stream is moving with velocity vj {zero order) then its susceptibility is 

Xj — — — 2 ■ ^ same direction) (5.169) 

{uj - kvj) 

Proof from equation of motion: 

-^E = ^ + v.Vv = {-iu + ik.Vj) V = -i (a; - kvj) v . (5.170) 

Current density 

j = PjVj +pj.v + pvj-. (5.171) 

Substitute in 

V.j + ^ = ik.vpj + ik.vp -iup^O (5.172) 

Pj = Pj , . (5.173) 



u — k.Vj 



Hence substituting for v in terms of E: 



X,eoV.E = p, = (5.174) 
rrij —I [uj — k.v,) 



which shows the longitudinal susceptibility is 



2 



X. = -^^-^ = ^^^ (5.175) 
TUjeQ (a;2 — kvj) {u — kvj) 

Proof by transforming frame of reference: 

Consider Galileean transformation to a frame moving with the stream at velocity v^. 

x^x' + Vjt ; t'^t (5.176) 
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exp i (k.x — Lot) — exp i (k.x' — (a; — k.Vj) t') 

So in frame of the stream, u' — u; — k.Vj. 
Substitute in stationary cold plasma expression: 



(5.177) 



2 2 

^ _ ^Pj ^Pj 



(5.178) 



Thus for n streams we have 



PJ 



[uj - kvj) 



2 ■ 



Longitudinal wave dispersion relation is 



e = 0. 



Two streams 



= e= 1 



UJ, 



UJ 



P2 



{uj — kvif {uj — kv2f 



(5.179) 

(5.180) 
(5.181) 



For given real k this is a quartic in a;. It has the form: 




Figure 5.12: Two-stream stability analysis. 

If e crosses zero between the wells, then 3 4 real solutions for u. (Case B). 
If not, then 2 of the solutions are complex: uj — ujr^ iuJi (Case A). 
The time dependence of these complex roots is 



exp {—icut) = exp {—icUrt ± ujit) . 



(5.182) 
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The +ve sign is growing in time: instability. 

It is straightforward to show that Case A occurs if 



\k{v2 -Vi)\ < UJ^i+UJ^2 



2 2 



3 
2 



(5.183) 



Small enough k (long enough wavelength) is always unstable. 

Simple interpretation {u!p2 -C ujpi, Vi — 0) a tenuous beam in a plasma sees a negative e if 

\kV2\ ^ UJpl. 

Negative e implies charge perturbation causes E that enhances itself: charge (spontaneous) 
bunching. 



Wave damping is due to wave-particle resonance. To treat this we need to keep track of the 
particle distribution in velocity space — > kinetic theory. 

5.9.1 Vlasov Equation 

Treat particles as moving in 6-D phase space x position, v velocity. At any instant a particle 
occupies a unique position in phase space (x, v). 

Consider an elemental volume rf^xrf'^v of phase space [dxdydzdvxdvydvz], at (x, v). Write 
down an equation that is conservation of particles for this volume 



5.9 Kinetic Theory of Plasma Waves 




= [v^f (x + dx±, v) - Va:f (x, v)] dydzd\ 

+ same for dy, dz 

+ [ttxf (x, V + dv^±) - ttxf (x, v)] d^^dvydvz 

+ same for dvy, dv^ 



(5.184) 




Figure 5.13: Difference in flow across a^-surfaces {+y + z). 



a is "velocity space motion", i.e. acceleration. 
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Divide through by (Px(Pv and take hmit 

= V.(v/) + V„.(a/) (5.185) 

[Notation: Use |^ ^ V; 

Take this simple continuity equation in phase space and expand: 

^ + (V.v) / + (v.V) / + (V,.a) / + (a.V,) / = 0. (5.186) 
Recognize that V means here ^ etc. keeping v constant so that V.v = by definition. So 

Now we want to couple this equation with Maxwell's equations for the fields, and the Lorentz 
force 

a=— (E + vAB) (5.188) 
m 

Actually we don't want to use the E retaining all the local effects of individual particles. We 
want a smoothed out field. Ensemble averaged E. 

Evaluate 

V^.a = V^.— (E + vAB) = — V^. (vAB) (5.189) 
m m 

= — B. (V^Av) = 0. (5.190) 

m 

So RHS is zero. However in the use of smoothed out E we have ignored local effect of one 
particle on another due to the graininess. That is collisions. 

Boltzmann Equation: 

9f df df (df\ 

\ / collisions 

Vlasov Equation = Boltzman Eq without collisions. For electromagnetic forces: 

|+v.?/ + ^(E + vAB)|^=0. (5.192) 
ot OX m av 



Interpretation: 

Distribution function is constant along particle orbit in phase space: ^/ = 0. 



dt^~ dt + dfdx^ dfdv 
Coupled to Vlasov equation for each particle species we have Maxwell's equations. 
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Vlasov-Maxwell Equations 

^ + v.|^^ + ^(E + vAB).|^^ = (5.194) 
at ax rrij ovj 

VAE = — , VAB = /xoj + ^^ (5.195) 

V.E = - , V.B = (5.196) 
Co 

Coupling is completed via charge & current densities. 

p = T.m-Y.^j j hd\ (5.197) 

3 J 

j j '' 

Describe phenomena in which collisions are not important, keeping track of the (statistically 
averaged) particle distribution function. 

Plasma waves are the most important phenomena covered by the Vlasov-Maxwell equations. 
6-dimensional, nonlinear, time-dependent, integral-differential equations! 

5.9.2 Linearized Wave Solution of Vlasov Equation 

Unmagnetized Plasma 

Linearize the Vlasov Eq by supposing 

/ = /o(v) + /i(?;)cxp2(k.x-cut) , /i small (5.199) 
also E = El exp i (k.x - cut) B = Bi exp i (k.x - cut) (5.200) 

Zeroth order /o equation satisfied by ^ — 0. First order: 

- zcu/i + v.zk/i + ^ (El + V A Bi) . ^ = 0. (5.201) 
m ov 

[Note V is not per se of any order, it is an independent variable.] 

Solution: 

/i = \ (El + V A Bi) (5.202) 

For convenience, assume /o is isotropic. Then ^ is in direction v so v A Bi.^ = 

— k.v) 

We want to calculate the conductivity cr. Do this by simply integrating: 

j = / g/ivd^t; = — / ^d?v .El. (5.204) 

J im J uj — k.v 
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= . ^ (5.203) 



Here the electric field has been taken outside the I'-integral but its dot product is with 
dfo/dv. Hence we have the tensor conductivity, 



g2 , ^,dfo 



[ ^ d\ (5.205) 
J uj — k.v 



im 

Focus on zz component: 

1 + X.. = 6.. = 1 + ^^ = 1 + ^ / -^^d\ (5.206) 
—lueo umeo J uj — k.v 

Such an expression applies for the conductivity (susceptibility) of each species, if more than 
one needs to be considered. 

It looks as if we are there! Just do the integral! 

Now the problem becomes evident. The integrand has a zero in the denominator. At least 
we can do 2 of 3 integrals by defining the 1-dimensional distribution function 

fziv,) = I f{w)dv,dvy (k = fcz) (5.207) 



r V ^ 

/ '^y- dv, (5.208) 
J UJ — kvz 



Then 

q2 . .,91. 

X ~ 

(drop the z suffix from now on. 1-d problem). 

How do we integrate through the pole at = |? Contribution of resonant particles. Crucial 
to get right. 

Path of velocity integration 

First, realize that the solution we have found is not complete. In fact a more general solution 
can be constructed by adding any solution of 

f+.f^O (5.209) 

[We are deahng with 1-d Vlasov equation: |£ + + = 0-] Solution of this is 

h^9{vt-z,v) (5.210) 
where g is an arbitrary function of its arguments. Hence general solution is 

<i j7"9/o 

/i = — '^&i^%{kz-ujt)^g{vt- z,v) (5.211) 
%\uj — kv) 

and g must be determined by initial conditions. In general, if we start up the wave suddenly 
there will be a transient that makes g non-zero. 
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So instead we consider a case of complex uo (real k for simplicity) where uj = uj^ + iuJi and 

LUi > 0. 



This case corresponds to a growing wave: 

exp{—iut) = exp(— ztUji + Uit) 



(5.212) 



Then we can take our initial condition to be /i = at t ^ — oo. This is satisfied by taking 
9 = 0. 

For uOi > the complementary function, g, is zero. 

Physically this can be thought of as treating a case where there is a very gradual, smooth 
start up, so that no transients are generated. 

Thus if UOi > 0, the solution is simply the velocity integral, taken along the real axis, with 
no additional terms. For 



Ui > 0, 



X 



ai 

D dv 



-dv 



umeo Jc uj — kv 

where there is now no difficulty about the integration because u is complex. 



(5.213) 



Figure 5.14: Contour of integration in complex f -plane. 
The pole of the integrand is at f = ^ which is above the real axis. 

The question then arises as to how to do the calculation if Ui < 0. The answer is by "analytic 
continuation", regarding all quantities as complex. 

"Analytic Continuation" of x is accomplished by allowing uj/k to move (e.g. changing the uJi) 
but never allowing any poles to cross the integration contour, as things change continuously. 



Remember (Fig 5.15) 



^ Fdz = ^ residues x 27ii 



(5.214) 



(Cauchy's theorem) 

Where residues = lim^.,^,, [F{z) / (z — Zk)] at the poles, Zk, of F{z). We can deform the contour 



how we like, provided no poles cross it. Hence contour (Fig 5.16) 
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Figure 5.15: Cauchy's theorem. 



Figure 5.16: Landau Contour 



We conclude that the integration contour for a;^ < is not just along the real v axis. It 
includes the pole also. 

To express our answer in a universal way we use the notation of "Principal Value" of a 
singular integral defined as the average of paths above and below 



P 



J V — vq 2 Uci Jc2\ 



-dv 



V - Vq 



(5.215) 



4 




Figure 5.17: Two halves of principal value contour. 



Then 



1' J 



V 



— ^dv 
u) — kv 



2 



(5.216) 



Second term is half the normal residue term; so it is half of the integral round the pole. 
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Figure 5.18: Contour equivalence. 

Our expression is only short-hand for the (Landau) prescription: 

"Integrate below the pole". (Nautilus). 

Contribution from the pole can be considered to arise from the complementary function 
g{vt — z,v). If g is to be proportional to exp{ikz), then it must be of the form g — exp[ik{z — 
vt)]h{v) where h{v) is an arbitrary function. To get the result previously calculated, the value 
of h{v) must be (for real cu) 



h{v) 



TT 



^1 dh 
mk dv 



5[v- 



(5.217) 



(so that 



-vgdv 



^ dfo 
k"^ dv 



■) 



(5.218) 



This Dirac delta function says that the complementary function is limited to particles with 
"exactly" the wave phase speed ^. It is the resonant behaviour of these particles and the 
imaginary term they contribute to x that is responsible for wave damping. 

We shall see in a moment, that the standard case will be uji < 0, so the opposite of the 
prescription Ui > that makes g = 0. Therefore there will generally be a complementary 
function, non-zero, describing resonant effects. We don't have to calculate it explicitly 
because the Landau prescription takes care of it. 



5.9.3 Landau's original approach. (1946) 

Corrected Vlasov's assumption that the correct result was just the principal value of the inte- 
gral. Landau recognized the importance of initial conditions and so used Laplace Transform 
approach to the problem 

A{p) = / e-'''A{t)dt (5.219) 







The Laplace Transform inversion formula is 

1 ps+ioo 

A{t) = — / ef"A{p)dp (5.220) 

where the path of integration must be chosen to the right of any poles of A{p) (i.e. s large 
enough). Such a prescription seems reasonable. If we make 3?(p) large enough then the A{p) 
integral will presumably exist. The inversion formula can also be proved rigorously so that 
gives confidence that this is the right approach. 
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If we identify p — > —iu, then the transform \s A — J e^*A{t)dt, which can be identified as 

the Fourier transform that would give component A oc e"*"^*, the wave we are discussing. 
Making 3?(p) positive enough to be to the right of aU poles is then equivalent to making '^{uj) 
positive enough so that the path in a;-space is above all poles, in particular oji > Q{kv). For 
real velocity, v, this is precisely the condition uui > 0, we adopted before to justify putting 
the complementary function zero. 

Either approach gives the same prescription. It is all bound up with satisfying causality. 



5.9.4 Solution of Dispersion Relation 

We have the dielectric tensor 



ujmeQ 



° -dv — m 



u — kv 



k'^ dv 



(5.221) 



for a general isotropic distribution. We also know that the dispersion relation is 

= (-A^2 + et)% = (5.222) 



-N^ + et 
-N^ + et 
e 



Giving transverse waves = et and 

longitudinal waves e = 0. 

Need to do the integral and hence get e. 

Presumably, if we have done this right, we ought to be able to get back the cold-plasma 
result as an approximation in the appropriate limits, plus some corrections. We previously 
argued that cold-plasma is valid if | ^ Vj. So regard ^ as a small quantity and expand: 



P 



Mi 

dr 



dv 



1 



dv 



1 + 



kv 



UJ 



kv 

UJ 



UJ J 



UJ 



n + 



2kv fkv' 

1 + + 3 — 

UJ \uj ^ 

3nT e 

+ ... 



+ 



dv 



dv (by parts) 



m uj^ 



(5.223) 



Here we have assumed we are in the particles' average rest frame (no bulk velocity) so that 
/ fovdv = and also we have used the temperature definition 

nT^ J mv'^fodv , (5.224) 

appropriate to one degree of freedom (1-d problem). Ignoring the higher order terms we get: 




re .uj'i dfo 

3 2+^^72^ ^ 

muj'^ k-^ n ov 



(5.225) 
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This is just what we expected. Cold plasma value was e = 1 — ^. We have two corrections 

1. To real part of e, correction = 3 (^^^ due to finite temperature. We could have 
got this from a fluid treatment with pressure. 

2. Imaginary part — > antihermitian part of e — > dissipation. 

Solve the dispersion relation for longitudinal waves e = (again assuming k real cu complex). 
Assume cui <^ cUr then 

,2 2 2 r T k"^ UJ^ldfo, , 

' mu^ k^ n dv 

- ^p{l + 3--2+-^-i:iT} (5-226) 



r 



1 2 .^rl^/0| 27ra;^l5/o| /croo'7^ 



For a Maxwellian distribution 

The difference between uo^. and uOp may not be important in the outside but ought to be 
retained inside the exponential since 



2rF 



1 + 3- 



mul 3 



uol I ( mujl 3\ 
So -.--,(«) #^-p(-^-2] (8-232) 

Imaginary part of uj is negative ^ damping. This is Landau Damping. 

Note that we have been treating a single species (electrons by implication) but if we need 
more than one we simply add to x- Solution is then more complex. 



5.9.5 Direct Calculation of Collisionless Particle Heating 

(Landau Damping without complex variables!) 

We show by a direct calculation that net energy is transferred to electrons. 
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Suppose there exists a longitudinal wave 

E = £; cos{kz - u;t)z (5.233) 

Equations of motion of a particle 

^ = —Ecosikz-ujt) (5.234) 

dt m ^ ' ^ ' 

Solve these assuming E is small by a perturbation expansion ^ = ^0 + ^1 + ..., z = ZQ{i) + 
Zi{t) + ... . 

Zeroth order: 

dvo 

= ^ = const , ZQ — Zi-\- vot (5.236) 

where Zi = const is the initial position. 
First Order 

E cos {kzo - out) ^ —Ecos {k {zi + Vot) - out) (5.237) 



dvi 


- s_ 


~dt 


m 


dzi 




~dt 


= Vi 



(5.238) 
Integrate: 

qE sin {kzi + kvo-ut) /r oQn^ 

Vi — ; + const. (5.239) 

m kvQ — ou 

take initial conditions to be Vi,V2 = 0. Then 

_ qE sin {kzj + Aut) - sin {kzj) (^oao) 
m Auj y ■ ) 



where Aa; = kvg — ou, is (-) the frequency at which the particle feels the wave field. 

qE 



Zi = 

m 



cos kzi — cos (kzi + Aout) sin kzi 
A^J2 ^ Aou 



(5.241) 



(using zi(0) =0). 

2nd Order (Needed to get energy right) 

dv2 qE [ , , . / , , . T 

— {cos {kzi + kvot — out + kzi) — cos {kzi + kvot — out) \ 
dt M 

qE 

= —kzi{- sin {kzi + Aojt)} {kzi<^l) (5.242) 
Now the gain in kinetic energy of the particle is 

]^mv^ - ]^rnvl = ]^m{{vo~^vi+V2 + ...f -vl] 

= ^{2,VqVi -\-vl-\- 2vqV2 + higher order} (5.243) 
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and the rate of increase of K.E. is 



d fl 2\ ( d.'^i d,Vi dv2 



We need to average this over space, i.e. over Zi. This wiU cancel any component that simply 
oscillates with Zi. 



^^0 



Vl 



dvi 
di 
dvi 
dt 



dV2 

dt 



(5.246) 



sin {kzi + Aa;i) — sin kzi 



cos {kzi + Aut) 



q^E"^ I sin {kzi + Aujt) — sin {kzi + Acut) cos Aujt + cos {kzi + Aujt) sin Aujt 



cos (kzi + Aut) ^ 

q'^E^ /sin Aa;t n . , \ 

^ cos^ (^2;^ + Aut) j 



\ Auj 
q^E^ 1 sin Acjt 

w? 2 Auj 



-q^E^ 



-kv, 







' cos kzi — cos ikzi + Aa;i) sin /c^j \ , , ^ , ^ 
^ ' ^ - t— — - 1 sin {kzi + Au;i) 



~q^ 

m? 
q^E^ kvo 



-kvn 



\ Acj2 
/sin Aa;t cos Acjt \ 



V Ac^2 
sin Aa;t 



Aa;2 



+ t- 



Au J 
cos Aa;f 



Aa; 



sin^ {kzi + Aujt) 



Auj 



(5.247) 



(5.248) 



Hence 



d 1 



-mv 



q^ 
2m 

q^E^ ^ 



sin Aujt 
Auj 



kvi 



sin Aujt 



2m 



-uj sin Aujt 



Ac^2 



+ 



Au;2 

ujt 



kvot 



cos Acjt 



Au; 



Auj 

cos Aa;i + i cos Aujt 



(5.249) 
(5.250) 



This is the space-averaged power into particles of a specific velocity vo- We need to integrate 
over the distribution function. A trick identify helps: 



—u , . ut ^ ^ 

- — - sm Aut + - — cos Aujt + t cos Aut 
Auj^ Auj 



d /a;sinAa;t 



dAu 



Auj 



1 J}_ 
k dvo 



UJ sin Aujt 
Auj 



+ sin Aujt 



+ sinAa;tj (5.251) 
(5.252) 
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Hence power per unit volume is 



d 1 



dt2 



mv^ ) f (vq) dvo 



2mk 
2mk 



d /a;sinAa;i 



u) sin Aa;t 



+ sin Aa;t J dvQ 

df 



+ sin Aojt I —^dvo 
J dvo 



(5.253) 



As t becomes large, sin Acut = sin(A;wo — uj)t becomes a rapidly oscillating function of vq. 
Hence second term of integrand contributes negligibly and the first term, 



CO sin Acut sin Acut 



oc 



Auj 



Aut 



(5.254) 



becomes a highly localized, delta- function-like quantity. That enables the rest of the inte- 
grand to be evaluated just where Aa; = (i.e. kvQ — uj — Qi). 




Figure 5.19: Localized integrand function. 



So: 



X = Aut 
and / '-^dz 



{kvo — u)t. 

7T SO 



fE^u^dl 
2mk k dv 



I 



smx 

X 



dx 



P^-E 



Trq^uj dfo , 



(5.255) 



(5.256) 



2mk'^ dv ^ 

We have shown that there is a net transfer of energy to particles at the resonant velocity 
from the wave. (Positive if | is negative.) 



5.9.6 Physical Picture 

Auj is the frequency in the particles' (unperturbed) frame of reference, or equivalently it is 
kvQ where v'q is particle speed in wave frame of reference. The latter is easier to deal with. 
Aujt — kv'ot is the phase the particle travels in time t. We found that the energy gain was of 
the form 

sin Aujt , , . s ,^ 

-d {Acut) . (5.257) 



/ 



Aujt 
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Cc-t ^1 




Figure 5.20: Phase distance traveled in time t. 

This integrand becomes small (and oscillatory) for Aut ^ 1. Physically, this means that 
if particle moves through many wavelengths its energy gain is small. Dominant contribution 
is from Au;t < n. These are particles that move through less than | wavelength during the 
period under consideration. These are the resonant particles. 




Figure 5.21: Dominant contribution 
Particles moving slightly faster than wave are slowed down. This is a second-order effect. 




Figure 5.22: Particles moving slightly faster than the wave. 

Some particles of this vq group are being accelerated (A) some slowed (B). Because A's are 
then going faster, they spend less time in the 'down' region. B's are slowed; they spend more 
time in up region. Net effect: tendency for particle to move its speed toward that of wave. 

Particles moving slightly slower than wave are speeded up. (Same argument). But this is 
only true for particles that have "caught the wave" . 

Summary: Resonant particles' velocity is drawn toward the wave phase velocity. 

Is there net energy when we average both slower and faster particles? Depends which type 
has most. 

Our Complex variables wave treatment and our direct particle energy calculation give con- 
sistent answers. To show this we need to show energy conservation. Energy density of 
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Figure 5.23: Damping or growth depends on distribution slope 



wave: 



W 



1 . 1 i^o, 1 

<sin^> Electrostatic Particle Kinetic 



+ n-m\v'^\ 
2 ' ' 



(5.258) 



Magnetic wave energy zero (negligible) for a longitudinal wave. We showed in Cold Plasma 
treatment that the velocity due to the wave is i; = Hence 

,2" 

(again electrons only) (5.259) 



W 



2 2 



2uM 



(5.260) 



When the wave is damped, it has imaginary part of a;, uji and 

dW — 1 dE"^ 

= W 

dt ^2 dt 

Conservation of energy requires that this equal minus the particle energy gain rate, P. Hence 

2 TT a; 1 (9/o 2 



_p I Wuj dfo I ^ 
2mk^ dv I fe 

2W^ ~ ^ 



1 



1 + 



UJ, 



X 



1 + 



(5.261) 



So for waves such that u ^ Up, which is the dispersion relation to lowest order, we get 

^TTUr 1 dfo 



UJi — UJ. 



P2k^n dv 



(5.262) 



This exactly agrees with the damping calculated from the complex dispersion relation using 
the Vlasov equation. 

This is the Landau damping calculation for longitudinal waves in a (magnetic) field-free 
plasma. Strictly, just for electron plasma waves. 

How does this apply to the general magnetized plasma case with multiple species? 

Doing a complete evaluation of the dielectric tensor using kinetic theory is feasible but very 
heavy algebra. Our direct intuitive calculation gives the correct answer more directly. 
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5.9.7 Damping Mechanisms 

Cold plasma dielectric tensor is Hermitian. [Complex conjugate*, transpose-^ = original 
matrix.] This means no damping (dissipation). 

The proof of this fact is simple but instructive. Rate of doing work on plasma per unit 
volume is P = E.j. However we need to observe notation. 

Notation is that E(k, uj) is amplitude of wave which is really 9ft(E(k, u) expi(k.x — ut)) and 
similarly for j. Whenever products are taken: must take real part first. So 



9ft (Eexpi (k.x - ut)) .3ft (jexpi (k.x - ut)) 



Ee'^ + E*e 



-i(f> 



k.x — iot.) 



E.je^*"^ + E.j* + E*.j + E*.j*e"^''^ 



(5.263) 



The terms e'^''^' & e'^^-?^ are rapidly varying. We usually average over at least a period. These 
average to zero. Hence 

(P) = ^[E.r + E*.j] = ^3ft(E.r) (5.264) 



Now recognize that j = cr.E and substitute 

1 



4 



E.<T*.E* + E*.<T.E 



But for arbitrary matrices and vectors: 

A.M.B = B.M^.A; 
(in our dyadic notation we don't explicitly indicate transposes of vectors). So 



E.<T*.E* = E*.cr*^.E 

1 



hence 



-E*. 



<j + cr 



.E 



(5.265) 



(5.266) 



(5.267) 



(5.268) 



If e = 1 H }^ — cr is hermitian e*"^ = e, then the conductivity tensor is antihermitian 



cr (if UJ is real). In that case, equation 5.268 shows that < P >= 0. No dissipation. 



Any dissipation of wave energy is associated with an antihermitian part of cr and hence e. 
Cold Plasma has none. 



Collisions introduce damping. Can be included in equation of motion 

m— = g (E + V A B) — mv u 
at 



(5.269) 



where p is the collision frequency. 
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Whole calculation can be followed through replacing m(— io;) with m(i/ — icu) everywhere. 
This introduces complex quantity in S, D, P. 

We shall not bother with this because in fusion plasmas collisional damping is usually neg- 
ligible. See this physically by saying that transit time of a wave is 



Size 



1 meter 



~ 3 X 10 seconds. 



Speed 3 x 10+^m/s 
(Collision frequency)"^ ~ lO/is — > 1ms, depending on Te,n, 



(5.270) 



When is the conductivity tensor Antihermitian? 



Cold Plasma: 



S 
iD 




-iD 
S 







p 



where 



S = 
D 



E 



PJ 

, ,2 



^ = 1 - E, 



PJ 
'J up- 



(5.271) 



This is manifestly Hermitian ij uj is real, and then a is anti-Hermitian. 

This observation is sufficient to show that if the plasma is driven with a steady wave, there 
is no damping, and k does not acquire a complex part. 

Two stream Instability 



(5.272) 



In this case, the relevant component is Hermitian (i.e. real) if both cu and k are real. 

But that just begs the question: If cu and k are real, then there's no damping by definition. 

So we can't necessarily detect damping or growth just by inspecting the dieletric tensor form 
when it depends on both uj and k. 

Electrostatic Waves in general have e = which is Hermitian. So really it is not enough to 
deal with e or x- We need to deal with a = —iuegX, which indeed has a Hermitian component 
for the two-stream instability (even though % is Hermitian) because cu is complex. 



5.9.8 Ion Acoustic Waves and Landau Damping 

We previously derived ion acoustic waves based on fluid treatment giving 



= 1 



CO, 



pe 



pi 



mini 



(5.273) 



Leading to o;^ ~ k'^ 



mi 
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Kinetic treatment adds the extra ingredient of Landau Damping. Vlasov plasma, unmagne- 
tized: 



k'^ Jc V — f dv n 

k 



dfoedv _ UJ^ 
fc2 



dfoi dv 



c V — dv n 



(5.274) 



Both electron and ion damping need to be considered as possibly important. 

Based on our fluid treatment we know these waves will have small phase velocity relative to 
electron thermal speed. Also Cs is somewhat larger than the ion thermal speed. 









\ ^ 

1 ^ 








\k - 







Figure 5.24: Distribution functions of ions and electrons near the sound wave speed. 
So we adopt approximations 



UJ . .UJ 



(5.275) 



and expand in opposite ways. 
Ions are in the standard limit, so 



X. £_ 
I — 9 



m Lo 



(5.276) 



Electrons: we regard r as small and write 



1 dv 

V — dv n 



- P 



1 dfoe dv 



V dv n 

2 rdfoe, 

dv 



n J 



n J dv"^ 
2 r m. 



2T, 



foedv for Maxwellian. 



me 



(5.277) 



Write Fo = fo/n. 

Contribution from the pole is as usual so 

,2 



pe 



rrie , . dFoe 



u)/k 



(5.278) 



139 



Collecting real and imaginary parts (at real cu) 



EricUr) — 1 + 



pi 



cot 



1 + 



3Ti 



m tot 



ei{ur) 



— TT 



2 9Foe 



pe 



dv 



pi 



dv 



\u!/k 



(5.279) 
(5.280) 



The real part is essentially the same as before. The extra Bohm Gross term in ions appeared 
previously in the denominator as 



^pi ^ ^ 

2 _ k'^Pili /,i2 



1 + 



3T,- fc2 



(5.281) 



Since our kinetic form is based on a rather inaccurate Taylor expansion, it is not clear that 
it is a better approx. We are probably better off using 



^ 

1 - 



1 



(5.282) 



Then the solution of £r(<^r) = is 



fc2 



rrii 



1 + k^X^r 



(5.283) 



De 



as before, but we've proved that 7e = 1 is the correct choice, and kept the k'^Xj^^ term (1st 
term of e^)- 

The imaginary part of e gives damping. 



General way to solve for damping when small 

We want to solve £(k, cu) — with cu — cUr + iuJi , cui small. 
Taylor expand e about real cUr'- 

de 

e{uj) ~ e{ujr) + iuJi-^\ujr 

aujr 

Let ujr be the solution of er{oJr) — 0; then 



d 

e{u) — iei{u!r) + icui- — £(0;^). 



This is equal to zero when 



Si M 

de{u!r) 



(5.284) 
(5.285) 

(5.286) 
(5.287) 
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If, by presumption, Ei <^ Sr, or more precisely (in the vicinity of £ = 0), dsi/dujr <^ 
then this can be written to lowest order: 



doJr 



Apply to ion acoustic waves: 



der (cOr) 
dujr 



a;. 



2 r 

pi 

r L 



2 + 4 



m,; UJ. 



so 



71 00'^ 



2 + 4^4 



pe 



dv 



o/k + ^. 



pi 



dv 



\Lo/k 



For Maxwellian distributions, using our previous value for a;^, 

" \2tiTJ 



dFoe, 
dv 



2 

e 2Te 



Te + m 



rrii 



rrir-. 1 + 



37i 



cxp 



me 


1 H 


3T/ 


2 




/l- 





V27r Vmj 

where the exponent is of order me/rrii here, and so the exponential is 1. And 



dv 



1 rrii 



1 + 



exp 



1 



371 



2r, 1 + A;2Ai, 



Hence 



UJr 




Te 1 + 



2T, 1 + k-'Xl 



T. 1 



1 m 



exp 



Te 



271 1 + k'^X^- 



De 



electron 



ion damping 



dSr/ duJr 

(5.288) 



(5.289) 



(5.290) 



(5.291) 



(5.292) 



(5.293) 



(5.294) 
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[Note: the coefficient on the first hne of equation 5.294 for ijJijijJr reduces to 
Ti/Te < 1 and kX^e < 1.] 

Electron Landau damping of ion acoustic waves is rather smalh ^ 
Ion Landau damping is large^ ~ 1 unless the term in the exponent is large. That is 



7r/8 for 



/ va^ 1 
70' 



unless 



> 1 



(5.295) 



Physics is that large ^ pulls the phase velocity of the wave: y = above the ion 
thermal velocity vu = \J~^. - If Cg ^ Vu there are few resonant ions to damp the wave. 

[Note. Many texts drop terms of order ^ early in the treatment, but that is not really 
accurate. We have kept the first order, giving extra coefficient 



1 + 



3T,; 



Te + 3T, 
Te + 6T, 



2T, 



(5.296) 



and an extra factor 1 + ^ in the exponent. When Tj ~ Tg we ought really to use full 

J- e 

solutions based on the Plasma Dispersion Function.] 

5.9.9 Alternative expressions of Dielectric Tensor Elements 

This subsection gives some useful algebraic relationships that enable one to transform to 
different expressions sometimes encountered. 



ujmeo Jc uj — kv w^meo k Jc \uj — kv J ov 



Q 



2 1 



_J_dfo 
meo k'^ Jc T — V dv 



dv 



1 Idfo 



k"^ Jc f — V n dv 

k 



dv 



UJ. 



p 
k^ 



P 



1 dFo^ .dFo, 

dv — vrz— — U 

— V 2v ov 



(5.297) 
(5.298) 

(5.299) 

(5.300) 



where = ^ is the normalized distribution function. Other elements of x involve integrals 
of the form 



dvi 



-d^V . 



(5.301) 



J UJ — k.v 

When k is in z-direction, k.v = k^Vz- (Multi dimensional distribution /q). 

If (e.g., Xxy) I 7^ z and j ^ I then the integral over vi yields / ^dvi = 0. If j = / 7^ 2; then 

dfo 



(5.302) 
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by parts. So, recalling the definition fz = J fdv^cdvy, 

Jo 



Xxx Xyy 



ujmto J — k.v 



< [ 

UJ J 



- — dv^. 
UJ J UJ — k.v 



(5.303) 



The fourth type of element is 



ujme. 



J UJ — kzV 



dh . 



(5.304) 



This is not zero unless fg is isotropic (= fo{v)). 
If / is isotropic 



dfo ^ dfodv^ ^ Vzdfo 
dvz dv dvz V dv 



(5.305) 



Then 



J UJ — 



dfo 
KzVz 



d^V 



f VxVz i 

J UJ — kzVz V 

/ — 1 — ^^"^ = 

J UJ — KzVz OVx 



(5.306) 



(since the Vj;-integral of dfo/dv^ is zero). Hence for isotropic Fo — fo/n, with k in the 
z-direction, 



X 



dvz 



C u>—kvz ^ 








— f Fo. dy 



0+ 
0+ 



^ k oj—kv, dvz 



(5.307) 



(and the terms 0"^ are the ones that need isotropy to make them zero). 



et 





(5.308) 



where 



UJ Jc UJ — kv 



-dv. 



u^l 



1 dFo, 



k^ Jc V — T dv. 



dv. 



(5.309) 
(5.310) 



k ^""^ 

All integrals are along the Landau contour, passing below the pole. 
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5.9.10 Electromagnetic Waves in unmagnetized Vlasov Plasma 

For transverse waves the dispersion relation is 

— = A^' = et = 1 - ^- / / , \ 5.311 

UJ-^ UJ n JC [OJ — KzVz) 

This has, in principle, a contribution from the pole at a; — kvz — 0. However, for a non- 
relativistic plasma, thermal velocity is ^ c and the EM wave has phase velocity ~ c. Con- 
sequently, for all velocities Vz for which foz is non-zero kVz -C uj. We have seen with the cold 
plasma treatment that the wave phase velocity is actually greater than c. Therefore a proper 
relativistic distribution function will have no particles at all in resonance with the wave. 

Therefore: 



1. The imaginary part of from the pole is negligible. And relativisitically zero. 
2. 



e* ~ 1 

= 1 

~ 1 

~ 1 



uj"^ n J-oo \ 



kvz k'^vl 
1 + + + .. 



dvz 



UJ 



e T 



2- + - 



1 + 



a;2 



(5.312) 



Thermal correction to the refractive index N is small because — ^ <^ 1. 



Electromagnetic waves arc hardly affected by Kinetic Theory treatment in unmagnetized 
plasma. Cold Plasma treatment is generally good enough. 
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